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(54) Adaptive intraf ield reducing of gaussian noise by fuzzy logic processing 



(57) The level of Gaussian noise in a memory field 
being scanned by rows is reduced by reconstructing 
each pixel by fuzzy logic processors, the latter process- 
ing the the values of pixels neighbouring the pixel being 
processed and belonging to a processing window 
defined by the last scanned row and the row being 



scanned, thus minimizing the memory requisite of the 
filtering system to a single row. The system perform an 
adaptive filtering within the current field itself and does 
not produce "edge-smoothing" effects as in prior adap- 
tive filtering systems operating on consecutive fields. 
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Description 

FIELD OF APPLICATION OF THE INVENTION 

5 The present invention relates to an adaptive filtering system to reduce the level of Gaussian noise in a memory field 

being scanned by rows. More in particular the system of the invention is adapted to reduce the noise of an image to be 
displayed on screen by implementing an adaptivity in response to brightness changes, according to a predefined sen- 
sorial response system. 

io BACKGROUND OF THE INVENTION 

Noise is inherent to any communication system and, in television broadcasting networks, originates in different 
ways. 

At the source of image, noise may be either generated by photon or thermal noise in a videocamera or by granular 
15 noise in films of cameras. Signal processing circuits and recorder play back systems add noise because the signal cir- 
culates in an intrinsically noisy loop. 

Finally, the signal may be carried to the transmitter via a distribution network that adds more noise, and broadcast 
via an intrinsically noisy transmission medium. The receiver circuits themselves add noise. Many techniques have been 
developed to reduce the problem of Gaussian noise and these are described in literature. The adaptivity of noise level 
20 to a certain image/video sequence is a key feature of the so-called "clever" filters. 

Generally, this adaptivity is obtained using non linear techniques such as averaging filters or filters with variable 
behavior depending on the estimation of some parameters. 

For instance, many known filters are obtained using scale and position estimates, or by estimating the temporal cor- 
relation (Recursive Averaging Filters) which is the most common techniques presently used in TV environments. 
25 These known techniques have intrinsic disadvantages and/or limitations. In the case of a video sequence strongly 
affected by noise, the filter may confuse the noise with motion effects resulting in a poor noise reduction effect. To better 
this situation, the user could automatically select amongst a set of preordered response characteristics such as to 
dynamically accentuate or reduce the filtering effect. However, these systems of arbitrary selection influences the tem- 
poral frequency, introducing a sort of "comet effect" which may becomes more annoying than an excessive edge 
30 smoothing. On the other hand, motion detection requires at least one field memory in the recursive loop of the filter with 
the consequent increment of total memory requisite of the system. 

SCOPE AND SUMMARY OF THE INVENTION 

35 It is evident the need and/or utility for a Gaussian noise filtering system that minimizes the memory requisite for 
noise filtering and that introduces an adaptivity to the local picture characteristics based on the sensorial response to 
brightness changes of the so-called Human Visual System without penalizing edge definition by an undue smoothing. 
The filtering system of the invention provides for an effective response to these important requisites by processing data 
relative exclusively to a certain memory field, thus minimizing the memory requisite. The essential characteristic of the 

40 system of the invention is the utilization of a fuzzy logic process which allows for a current definition of the adaptivity 
parameters bearing in mind the characteristics of a certain scheme of sensorial response, such as for example the 
scheme of response to brightness changes known as the Human Visual System (HVS). The method of the invention 
permits to curt the noise affecting a certain area of a field scanned by rows, without employing information devided from 
a preceding field or picture which, by contrast, is the classic approach of the known systems wherein the noisy pixel is 

45 substituted with an averaged pixel value of the nearby ones belonging to a certain processing window. 

This known approach has indeed the serious disadvantage of smoothing the edges and the contours of the objects 
depressing image definition. 

By contrast, according to the present invention, each noisy pixel singled out within a certain processing window, is 
substituted with a weighted average of the surrounding pixels, according to a prestablished degree of similarity between 

so the same pixel to be substituted and each of the nearby ones. 

In this manner, if K ( provides a measure of the degree of similarity between the pixel X to be filtered and a generic 
pixel X; belonging to the process window, the output filtered pixel value is given by: 

X(') ou , = K i *X + (1-K i )*X i (1) 

55 

If K| is equal to 0, X and X, are equal and thereby X may be substituted with X ( . Contrarily, if K, is equal to 1. X and 
X| are radically different from each other and in this situation their exchange would introduce a serious error. Therefore, 
according to the present invention, the filter output is obtained through a superimposition of the effects of several ele- 
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mentary filters whose function is described by equation (1), and such that the output pixel is expressed by the following 



X represents the value relative to the gray level of the pixel being processed; 

Xj represents the value relative to the gray level of a pixel of neighbouring the pixel under processing 

Kj provides the measure of similarity between the pixel X and the adjacent pixel Xj, in practice K, represents the 
value produced by a fuzzy logic processing block according to a membership function corresponding to the sen- 
tence "the difference between the pixels is low". 

N is the number of neighbouring pixels considered. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 shows a processing window 

Figure 2 is an example of the membership function used to weigh the computed differences between the pixel 
under processing and the neighbouring pixels within the processing window of Figure 1. 

Figure 3 depicts a suitable membership function for generating the processing coefficients of a weighed average. 
Figure 4 is a basic scheme of a Gaussian noise filtering system according to the present invention. 
Figure 5 shows a block diagram of a circuit for defining a certain processing window. 

Figure 6 is a diagram of the computation blocks of the absolute values of the differences and extraction of the min- 
imum difference and of the maximum difference. 

Figure 7 is a functional diagram of the fuzzy logic processing block of the pixel being processed and of the minimum 
and maximum differences to generate a pair of adaptive threshold values. 

Figure 8 is a functional diagram of the plurality of fuzzy logic processing blocks of each computed difference to gen- 
erate a set of weighed coefficients that affect an identical number of elementary filters and of the adder block regen- 
erating the output pixel value. 

DESCRIPTION OF THE METHOD OF THE INVENTION 

Let us consider a processing window as the one shown in Figure 1 , being X the pixel to be processed. The pixels 
X1, X2, X3 are the pixels of a scan line preceding the current scan line to which the pixel to be processed belongs neigh- 
bouring the pixel X, whereas X4, X5, X6 and X7 represent neighbouring pixels of the pixel X, belonging to the current 
scan line. 

According to a fundamental aspect of the invention the system must preliminarily carry out an estimation of the 
noise level. 

For this purpose, the system computes the absolute values of the differences D, between the pixel to be filtered. X 
and the neighbouring ones, X h and finds the maximum difference D max and the minimum difference D min among the 
absolute values of the computed differences Dj. 

Using a fuzzy logic processing, different membership functions can be found that fulfill the sentence: "the absolute 
difference value |D,| is smair. but if. by hypothesis, the produced value K, varies according to the following law: 



relationship: 



X out = X X(i) out = 1 +[ I K i ' X + ( 1-K,) * X,\ 
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several drawbacks would arise: 

i) if the processing window does not contain edges, or boundaries of objects or features, D max and D min will be 
exclusively correlated to the noise, which therefore would not be significantly reduced; 
5 ii) on the contrary, if the processing window contains edges or boundaries, these would be unduly smoothed 

because of a certain loss of edge information. 

To avert the above mentioned drawbacks, a membership function based on the estimation of the noise level is 
implemented. 

io If the filter processing window coincides with a relatively homogeneous area of a picture, all the differences D| will 
be close in terms of their absolute value and correlated to the noise. Thence, it is possible to select a membership func- 
tion that is able to generate threshold parameters Thl and Th2 that substantially coincide with D max . 

On the contrary, in case the pixel to be filtered belongs to an edge or boundary region of the image, the filtering 
action is made significant only where the differences among the pixels of the processing windows are small, that is to 
is say, within the interval of the differences |Dj| defined by [D min ;(D min +D max )/2], 

The threshold parameters Th1 and Th2 are then computed according to the following equations: 

Th1 =K n *D max + (1-K n )*D min (4) 

20 

T"^K n 'D max + (1-K n )- Dma ^ Dmi " (5) 



It is easy to demonstrate that Th1 is always greater or equal to Th2. 
25 The range of the fuzzy variable "Low" may then be defined in terms of D max , D min , and K n , where K n influences the 

filter behavior according to the noise level determined for the area of the image defined by the processing window. 

As far as K n is concerned, in case the pixel to be processed belongs to a relatively homogeneous region of a picture 
affected by noise. D max and D mjn will be mainly correlated to the noise content and will have not so different values. In 
other words, this means that the difference (D max -D min ) is relatively small. 
30 If the pixel X to be processed belongs to an edge or boundary line, the difference (D max -D min ) will be large. Hence, 
the value of K n , is obtained, according to the method of the invention, by using a weighed function of the type depicted 
in Figure 2. 

The parameter A depends on the brightness of the pixel to be processed. For example, it can be established on: 

35 A = x/A. (6a) 

A = (255-x)/A. (6b) 

By way of an example, a satisfactory value for X may be 4. 
40 in choosing the parameter A for equation (6a), a standard scheme of sensorial response to contrast as the so- 
called Human Visual System (HVS) may be taken into consideration. 

According to Weber's law, equal increments on the logarithm scale of pixel luminance value are perceived to be 
equally different among from one another therefore: 



wherein AC represents the contrast variation. 
This means that: 



A(log)ocAC (7) 



Ax/x = a (8) 



where a is a constant. 

Hence, when an edge region is encountered, the condition axxix is complied with. 

While Weber's law is generally used to detect the minimum gray level to be quantized without introducing any false 
55 contour effect, according to the method of this invention, Weber's law is used instead to modulate the filtering action in 
function of the luminance characteristics of the pixel being processed. If the parameter A is chosen according to the 
relationship (6b), a stronger filtering action is implemented in dark areas than in bright areas of the image, if compared 
to that would be obtained by using equation (6a), this in consideration of the fact that the human eye is more sensible 
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to variations due to Gaussian noise in dark areas than in bright areas. 

DESCRIPTION OF THE ARCHITECTURE OF THE INVENTION 

A functional block diagram of a system implementing the above described technique is illustrated in Figure 4. 

After "constructing" the required processing window, according to common buffering techniques, which obviously 
involve information relative to the pixels belonging to at least two successive scanning lines or rows, a dedicated DIF- 
FERENCE COMPUTATION circuit computes the differences D f (D1, D2 D7) between the pixel luminance values 

(X1 . X2 X7) of the processing window neighbouring the pixel X being processed. 

The block, D max & D M(N DETECTION, identifies and selects by making them available through two dedicated out- 
puts, the maximum and minimum values of the differences D, input to the block. 

A first fuzzy logic Adaptive Threshold Detection processing block receives as input the pair of values D max and D m j n 
of the current set of differences D| and computes the parameter K n according to the following function referred to as 
MEMBERSHIP 1 , shown in Figure 2 and expressed by: 

if x > A => K n = 0 

ifO<x<A=>1 - (x/A) 

generating a corresponding pair of threshold values Thl and Th2, according to equations (4) and (5), respectively. 

A second fuzzy logic LOCAL ADAPTIVE NOISE SMOOTHING processing block carries out the independent 
processing of each of the differences Dj, using the previously generated threshold parameters Th1 and Th2 according 
to the function referred to as MEMBERSHIP 2, which may be depicted as done in Figure 3, and defined by: 

if | D j | < Th1 => K j = 1 
if | D j | > Th2 => K j = 0 

if Th 1 * | D ,| * Th2 =, K , = J ^A Hl - TH I^ m - 2 

producing a set of output values K it each representing the measure of the degree of similarity between the pixel to be 
filtered, X, and a respective neighbouring pixel X, and inside the processing window, according to equation (3). 

A final filtering block separately processes each of the difference values Dj in function of the respective weighing 
coefficient K, and, through a final adder stage I, produces the X,-, ut pixel as the sum of the weighed contributions, 
according to equation (2). 

The processing window can be easily realized by a circuit as that shown in Figure 5, whose functioning is immedi- 
ately comprehended by a person skilled in the art. 

It is evident that the memory requisite is substantially limited to the information of a simple scan line or row (Line 
Memory) and therefore it is a minimum fraction if compared to the memory requisite needed to store a whole field or 
picture. 

A functional diagram of the block DIFFERENCE COMPUTATION and D MAX & D M(N DETECTION is illustrated in 
Figure (6). 

Even in this case the circuit functioning is readily comprehended by a skilled person. 
The ADAPTIVE THRESHOLD DETECTION block is schematically detailed in Figure 7. 

According to a preferred embodiment, a fuzzy logic processing block according to a certain membership function 
(MEMBERSHIP 1), receives as input the difference values D min and D max , the value of the pixel under processing and 
a certain parameter X and outputs a value K n which, applied to the threshold generating blocks TH1 Gen and TH2 Gen, 
conditions the generated threshold values Th1 and Th2. 

Structural details of the LOCAL ADAPTIVE NOISE SMOOTHING and FINAL FILTERING block are shown in the 
functional diagram of Figure 8. 

The plurality of fuzzy logic circuits (MEMBERSHIP 2), generate a set of weighing coefficients Kj, by processing 
each difference value according to the threshold pairs Th1 and Th2 previously generated by the circuit of Figure 7. The 
membership function used may be the one depicted in Figure 3. 

The computation of the set of coefficients Kj allows to curt the noise level through the FINAL FILTERING block by 
using a number of partial filters suitable to produce a corresponding set of weighed contributes X(i) out that are then 
added by the final adder z which outputs the correct pixel Xo Ut . 
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Claims 

1. A method of filtering Gaussian noise contained in a data stream relative to a certain field scanned by rows, char- 
acterized in that it comprises 

5 

storing the information of a preceding scanned row defining together with the data of a row being currently 
scanned a certain bidimensional processing window of a certain number of pixels; 

computing the difference (Dj) between the value of a pixel (X) of the row being scanned and the value of the 

10 neighbouring pixels (X1 , X2 X7) belonging to said processing window and determining the minimum (D min ) 

and the maximum (D max ) of the computed differences; 

processing by fuzzy logic the values of said pixel (X), of said minimum and maximum differences (D min , D max ) 
and of at least a parameter (X), derived from a certain sensorial response (HVS), according to a first member- 
's ship function (MEMBERSHIP 1), generating a first weighing coefficient (K n ); 

generating a pair of threshold values (Th1 , Th2) of said differences according to said first weighing coefficient 
(K n ); 

20 processing by fuzzy logic each of said computed differences (Dj) and said pair of threshold values (Th1 , Th2), 

according to a second membership function (MEMBERSHIP 2), producing a corresponding set of weighing 
coefficients (Kj); 



25 



30 



combining each of said computed differences (Dj) with a respective weighing coefficient (K|); 

adding the results of said combinations producing a corrected value (X oul ) of said pixel (X). 

2. The method according to claim 1 . characterized in that said first weighing coefficient (K n ) is generated by using a 
membership function (MEMBERSHIP 1) defined by: 

if x > AK n =0 

ifO <x<A=> K n = 1 -(x/A) 

35 3. The method according to claim 1, characterized in that said pair of threshold values (Th1, Th2) is generated 
according to the following equations 

Th1 = K n *D max + (1-K n )*D min 
D „„ + D . 



40 



Th2 = K n *D max + (1-K n )* 



max mm 



4. The method according to claim 1 , characterized in that said weighing coefficients (kj) are generated using a mem- 
bership function (MEMBERSHIP 2) defined by: 

if | Dj | < Th1 => K, = 1 
if | D, | > Th2 => Kj =0 

I D, I 7772 



if Th1 < | Dj | < Th2 => Kj 



Thl - 7772 7771 - 77)2 
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5. A system of adaptive filtering to curb Gaussian noise within a same field of an image, characterized in that it com- 
prises 



6 



EP 0 878 776 A1 



buffering circuits of respective data sequences relative to a certain row of said field being scanned and to at 
least a previously scanned row of the same field, capable of defining a bidimensional processing window of pix- 
els belonging to said consecutive rows (X. X1, X2, X3, X4, X5, X6, X7); 

a circuit computing the differences (Dj) between the value of a certain pixel (X) and the neighbouring pixels 
(X1 X7) of said processing window and that determining the maximum and minimum computed difference 

a first fuzzy logic processor (MEMBERSHIP 1 ) adaptively generating a pair of minimum and maximum thresh- 
old values (Th1 . Th2) of said computed differences (Di). depending on the value of said pixel (X), of said min- 
imum and maximum difference (D mjn , D max ) and of at least a parameter (X) derived from a certain predefined 
sensorial response (HVS). according to a first membership function; 

a plurality of second fury logic processors (MEMBERSHIP 2), each combining a respective difference (Dj) with 
said pair of threshold values (Th1 , Th2), according to a second membership function and generating a respec- 
tive weighing coefficient (kj); 

a plurality of combinating logic circuits (PARTIAL FILTER^) each modifying a respective difference value (Dj) 
according to the respective weighing coefficient (k,); 

an adder circuit adding the modified values of said computed differences (Dj) and outputting a corrected pixel 
value (X out ) of said pixel (X). 
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